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Abstract 

Photometric and spectroscopic observations and analysis 
of the eccentric eclipsing binary V459 Cassiopeiae 
(e = 0.0244) were performed by Lacy et al. (2004). 
Observations of minimum light show the presence of 
apsidal motion. In order to find the observed rate of 
apsidal motion, I followed the procedure described by 
Guinan & Maloney (1985). A new observed rate of 
apsidal motion of 15°. 2/100 yr with a period of 2400 
yr is computed which is not in agreement with the one 
reported earlier. Also the advance of the periastron 
is calculated theoretically by taking into account the 
Newtonian (classical) and general-relativistic effects 
according to the physical and orbital parameters of the 
system. The theoretical value of 2°. 64/100 yr is obtained 
which is 5.75 times smaller than the observed rate of the 
apsidal motion. 

Key words: stars: binaries: eclipsing: V459 
Cas - stars: apsidal motion 



1 Introduction 

V459 Cassiopeiae (BV5, GSC 04030-01001; V„,ax = 
+10.33; P=8.45 days) is a double-lined eclips- 
ing binary where its variablity was discovered by 
Strohmeier (1955). A measurement on the eccentric- 
ity of the system was done by Busch (1976) after Mei- 
nunger & Wenzel (1967) recognized the eccentricity 
of its orbit. Through a series of spectroscopic and 
photometric observations (in UBV and ubvy) car- 
ried out from 1984 to 2004, the absolute dimensions 
of the physical and orbital parameters of V459 Cas 
was obtained after precise, simultaneous analysis of 
the light and velocity curves of the system (Lacy et 
al., 2004; hereafter L04). According to L04, the cor- 
responding radii and masses of the components are 



Ri = 2.009 ± O.O13i?0, i?2 = 1-965 ± O.OlSi?© and 
Ml = 2.02 ± 0.03Mq, M2 = 1.96 ± 0.03Mq respec- 
tively with an eccentricity of e=0.0244. But the rate 
of apsidal motion, lo, or the period of apsidal rev- 
olution, J7, is not well determined which is an im- 
portant parameter for probing the stellar structure 
of the eclipsing binary components. I investigate the 
observed rate of apsidal motion in Sec. 2 by following 
the method of Guinan & Maloney (1985) (hereafter 
GM85) to measure the observed rate of apsidal mo- 
tion. The comparison between the theoretical 
and the observational values of w is presented in Sec. 
3. The final results and conclusions are given in Sec. 
4. 



2 Observed rate of apsidal mo- 
tion 

Due to the deep, narrow eclipses of V459 Cas the 
rate of apsidal motion can be determined by analysis 
of primary and secondary eclipse timings. In their 
paper, GM85 described the procedure that must be 
followed to determine the apsidal motion rate from 
the change in the displacement of the secondary min- 
imum from the half point (0.5 phase) according to 



D ^ {t2- h) - O.bxPeriod, 



(1) 



where ^2 and ti are times of secondary and primary 
minima, respectively and D is related to the longi- 
tude of periastron lo by the formula given by Sterne 
(1939,a,b) 
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In Eq. ((SJ, e is the orbital eccentricity and P is 
the orbital period in days. Fig. 1, presents the (0-C) 
diagram from all the observed times of primary and 
secondary minima. To plot the (0-C) diagram for 
the primary minima (upper panel of Fig. 1), I used 
the times of minima listed by L04 as well as two other 
minima reported by Lacy et al. (1998) and Hiibscher 
(2005). The residuals are calculated according to the 
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Figure 1: Differences between tfie observed and the 
computed (0-C) times of primary and secondary 
minima, computed from Eq. Q, versus Julian Day 
number and year. The observed times of secondary 
minima are hsted in Table 1. 



ephemeris given by L04 

Min.I = HJD 2452565.67170-1- 8.45825381d x E. 

(3) 

For all the secondary minima in Table 1 , D and from 
that the uj are calculated using Eq. (jSJ. As you can 
see from Table 1, the slow increase (or decrease) in D 
(or Lu) is due to the advance of the line of apsides of 
the orbit and by computing the slope of a line fitted to 
the secondary minima on lo versus Julian Day number 
diagram, I determined a new rate of apsidal motion 
^obs ^ ± 6.75 (see Fig. 2). 

It must be noted that among the observed sec- 
ondary minima listed in Table 1, the error values of 
the last five minima are at least 100 times smaller 
than the other minima. Therefore, I use only these 
five minima to measure the tj"''^. Also in some cases 
the error in secondary minima are so high and Eq. 
doesn't give any lu correspond to the D, calcu- 
lated from Eq. (QJ since as it is obvious from Fig. 
3, for any set of values of e and P in Eq. (0), -D is 
always limited between its maximum and minimum 



Figure 2: Plot of the longitude of periastron {oj) ver- 
sus Julian day number according to the displacement 
of observed secondary minima listed in Table 1 from 
a half period point (D). The inset shows only the 
last five minima of Table 1 for which the calculated 
slope of the diagram is based only on these five sec- 
ondary times of minima. A least square fit to the 
observations, yielding a slope cj°^* — 15.18deg/100 
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Figure 3: D as a function of lu using Eq. 
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Figure 4: Relative orbit of V459 Cas, drawn to scale, 
from the quantities given by Lacy et al.(2004). 



values (e.g. Dmi„ < D < D^ax)- Thus, I can not 
find any uj correspond to Z? > Dmax or D < Dmin- 
Finally Fig. 4 shows the relative orbit of V459 Cas, 
drawn to scale. Actually what can be measured from 
Eq. ^ \s UJ and according to Fig. 4, w"'"' < but 
L04 have considered uj' as the angle of line of apside. 
One can reveal from Fig. 4 that lo = 2ii — lj' and 
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3 Comparison with theory 

Basically the theoretical rate of apsidal motion is due 
to the contribution of two terms; a classical term as 
well as the general-relativistic term. Then the total 
rate of apsdial motion is equal to 



•theo 
^tot 
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GR ■> 



(4) 



where io*cl° denotes the classical or Newtonian 



term and Wq^' 



uJq^" is the relativistic contribution. I am 
going now to calculate each term in equation Q with 
the assumption that the rotation equators of both 
stars are coincident with their orbital plane. In classi- 
cal mechanics, the advance of the line of apsides arises 
due to the lack of spherical symmetry in the shape 
of the stars. Tidal and rotational distortion are two 



mechanisms lead to this asymmetry which depend 
on the fractional radii, the internal mass distribu- 
tion, and the axial rotation of the stars. ojq'I" can be 
determined using the expressions given by Cowling 
(1938), Sterne (1939), and Kopal (1959) 
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where /2(e) = (1 + \e^)(X - (?Y^ , P is the orbital 
period in days, M\ and Mi are the masses of stars 
in terms of M©, A;2,i and fc2,2 are known as the 
apsidal motion constants of the component 1 and 
2 respectively, r\ and ri are the fractional radii of 
the stars (r — ^\A= semimajor axis of the orbit), 
(Dr,i and ijJr,2 are the star's angular rotation speeds, 
and is the mean angular Keplerian velocity equal 
to In series expansion /2(e), I did not keep 

the terms involving and higher orders since they 
have not any significant contribution to the result 
of summation in /2(e). L04 indicated that the 
components of V459 Cas are main-sequence stars 
with an age of about 525 Myr. Therefore in order 
to calculate ki which is correspond to the second 
harmonic of the mutual tidal distortion, I used 
Table 1 of Jcffery (1984) where he has computed 
^2 from modern stellar-interior models of evolving 
main sequence stars. From the estimated values 
for chemical composition (X — Q.lVl^Z — 0.012), I 
adopt a nearly equal value of k^^i ~ ^2, 2=0. 0042 for 
both components. Using the profiles of absorption 
lines in the star's spectrum, it is possible to estimate 
the angular rotation velocities Cor of the stars from 
their projected rotational velocity, 14 sin i. L04 
have estimated Vr sin i from the absorption lines 
of Mg//A4481 in the spectrogram of V459 Cas. 
From these measurements, V^sini=54 kms~^ and 



43 kr 



for the first and the second compo- 



nents respectively. Therefore (ir.i/'i'fc=4.49 and 



3 



Wr,2/wfc=3.66 are obtained from LJr,i = Vrot/Ri- 
Substituting all of the parameters and the values 
from L04 into equation I find 



' .theo 



1°.22± 0.12/100 yr. 



The expression for the advance of the line of 
apsides due to general relativity has a simple form 
and according to Levi-Cevita (1937) and Kopal 
(1959) is equal to 



^hn'-^)- 9.2872 . 10-^^Ml±^ 



(Ml 



(6) 



where Mi and M2 are in Mq and P is in days. 
Putting again the value of the parameters into 
equation 10 yields 



theo 
GR 



1°.42± 0.01/100 yr. 



Thus the combined classical and general rela- 
tivistic rate of apsidal motion is equal to 



theo 
CL+GR 



2°.64± 0.12/100 yr, 



which is calculate according to equation 
To calculate the standard deviations, I used the 
errors assigned to each parameter given in the refer- 
ences. The above value indicates that in comparison 
is about 5.75 times smaller with a 
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discrepancy of 12°.54deg/100 yr. 

4 Results and discussion 

V459 Cas has been observed and analyzed carefully 
by L04. In this paper I followed the method of GM85 
to measure a value for the observed rate of apsidal 
motion. Table 2 presents the results of independent 
determinations of the rate of apsidal motion together 
with the corresponding period of apsidal revolution 
U . The relation between the apsidal motion period 
U and w"'"* has a simple form 



U 



360P 

'• ; 

^ohs 



(7) 



where Cuobs is expressed in degrees per cycle and P 
is the anomalistic period expressed in days. Also the 
observed rate of apsidal motion is 5.75 times greater 
from the expected theoretical value. It is not the 
first time that such a discrepancy is observed. For 
example in the case of VI 143 Cygni the observed 
rate of apsidal motion is about one-third of the value 
estimated from the theory but according to the new 
studies it seems that this deviation is decreasing with 
the expansion of observational time (Dariush et al., 
2004). On the other hand in the case of highly eccen- 
tric eclipsing binary, DI Herculis, the observed rate 
of apsidal motion is one-seventh of the estimated the- 
oretical value. Tough a series of investigations have 
been made during the past decades to find a clear 
explanation for this high discrepancy (GM85; Claret, 
1998; and Claret & Willems, 2002) but the case of 
DI Herculis is still unsolved. The method of GM85 
is expected to gives the value of exactly but to 
determine a more accurate value for we need 
more accurate timings of secondary minima since up 
to now, our observational baseline covers only a very 
small fraction of the period of apsidal revolution U . 
Therefore it seems that the reason for such a large 
discrepancy between the theory and the observation 
is mainly due to the small number of observational 
data points. 
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Table 1: The photoelectric times of secondary min- 
ima for V459 Cas, together with the values of D and 
w, computed using Eqs. and (|2Jl according to the 
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Table 2: Determined rate of apsidal motion for V459 
Cas. 



^obs (deg/100 yr) 




Source 


6.04 ±4.74 


6100 


L04 


15.18 ±6.75 


2400 
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